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Abstract 

A recursion operator is constructed for a new integrable system 
of coupled Korteweg-de Vries equations by the method of gauge- 
invariant description of zero-curvature representations. This second- 
order recursion operator is characterized by unusual structure of its 
nonlocal part. 



1 Introduction 

Recently, the Painleve classification of coupled Korteweg-de Vries (KdV) 
equations was made in [Q, and the following new system possessing the 
Painleve property was found there: 

_ _5ci 364 _M _M 

"^t TTT^xxx , , "^xxx C\UxU ; UxV ; VxUi 

12 Ml di di 

6di 76ci 6c? 

Vt = —Uxxx —Vxxx + diUxU - ClUxV - ZCiVxU - —VxV, (1) 

where 6, C\ and d\ are arbitrary nonzero constants. 

More recently, the coupled KdV equations (P appeared in ||^, also as a 
system possessing the Painleve property, in the following form: 

= Uxxx ~ ^Uxxx ~ ^'^UUx + ISvUx + ISuVx, 

Vt = Uxxx — 5Vxxx — I2uux + I2vux + 6uVx + ISvVx, (2) 



which is related to (|I|) by the transformation (|T])t-^(|2|): 

U ^ OU f , V f , t 1-^ - — t. [6] 

2 ' 2ci ' 5ci ^ ' 

It was found in that the system (|^) admits the zero- curvature represen- 
tation (ZCR) 

DtX = D,T-[X,T], (4) 
which is the compatibihty condition for the overdetermined hnear system 

= Xm, = T^, (5) 

with 

(cr u — v crdw — 4t>)\ 
I -2(T M - , (6) 
Of a I 

where a is an essential ('spectral') parameter, Dt and D^. stand for the total 
derivatives, the square brackets denote the commutator, \E' (x, t) is a three- 
component column, and an explicit (but cumbersome) expression for the 3x3 
matrix T can be seen in p| . 

In the present paper, we study this new integrable system of coupled KdV 
equations in its equivalent form, which is characterized by more transparent 
structure of X in its ZCR. Namely, making the transformation 



1 4 1 1 

2u f , V ^ -u f , t ^ — t, 

3 3 3 2 ' 



m^S-^, X^SXS-\ T^STS-\ ^ = diag 0,1,0 



(7) 



in (H), (H) and (H), we obtain the system 

Ut = 'iUxxx — Vxxx — 12uUx + VUx + 2uVx, 
Vt = 9Uxxx — '^Vxxx — l^VUx — 6uVx + 4f f^. 



which admits the ZCR (HI) with 




X = \1 -2a u (9) 
\0 1 a J 

(the explicit expression for T is omitted, since not used in what follows). 
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Our aim is to obtain a recursion operator for the system (^. We derive 
the recursion operator from the matrix X (^) of the system's ZCR. We do 
this, following the method of gauge-invariant description of ZCRs of evolu- 
tion equations, introduced in (note also a more general gauge-invariant 
formulation of ZCRs, developed in 0]). 

The paper is organized as follows. In Section 2, we find the complete class 
of coupled evolution equations admitting ZCRs (§) with the matrix X given 
by (y). This class is determined implicitly, through a three-term recurrence, 
which is explicitly solved in Section 3, giving rise to the recursion operator 
sought. In Section 4, we prove that this operator satisfies the sufficient 
condition of being a recursion operator in the sense of symmetries. Section 5 
contains discussion of the results. 



2 Class of systems 

Let us find the complete class of coupled evolution equations 

f (3^, t , U, . . . , Ux...xy • • • ; ^X...X^ ; 

Vt = gix,t,U,..., Ux...x, V,..., Vx...x) (10) 

which admit ZCRs (^) with the matrix X given by (^ and without any 
restrictions imposed on the 3x3 matrix T{x, t,u, . . . , Ux...x, v, . . . , Vx...x, cr). 
We solve this problem algorithmically, using the technique of cyclic bases of 
ZCRs, developed in ^ (see also ^ for more examples). 

Since the matrix X does not contain derivatives of u and v, the 
characteristic matrices of the ZCR are simply Cu = dX/du and = dX/dv: 
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a = 1 , a = . (11) 



Using the operator V^;, defined as S/xY = DxY — [X, Y] for any (here, 3x3) 
matrix Y, we compute VxC^, VxCy, V^.C„, V^Cy, etc., and find that the 
cyclic basis (i.e. the maximal sequence of linearly independent matrices V^.C„ 
and V^.C^, i,j = 0, 1,2, . . .) is 

{Cu, Cy, VxC'u, V^C*!,, V^.Cti, V^Ct,, V^^Cu, V^.Ctt} , (12) 
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and that the closure equations for this cychc basis are 

V^C^t = aiCu + CL2Cv + a-^S/xCu + 0'4:^xCv + asV^C^ + agV^C^, 

+ brVlCu + bsVtCu, (13) 

where 

0-2 = — 6Mxxx + ^Vxxx + 54:UUx — 3vUx — 22uVx + "ivvx 

45 

as = '^Uxx - 22m^ + ?,uv - 99ct^m + yO"^^ - 81a^, 

Q 

04 = — 24M2,.a. + IVxx + 54m — 36ui7 + -t> + 54cr n — 9o- 

05 = bux, ae = -36^2; + Sv^, 

3 2 
07 = 13n — -f + 18cT , ag = 0, 

oi = -Ux, 02 = -u^ + Vx, b^ = u + -a , 
3 1 

64 = -M+2^' ^T = -2' h = b6 = bs = 0. (14) 

Thus, in the case of the matrix X (^, the cychc basis is eight-dimensionaL 
Therefore the singular basis (see [§] for the definition) can be one-dimensional 
at most; in fact, it consists of the unit matrix / = diag (1, 1, 1) and, due to 
V^/ = 0, has no effect on the form of the represented equations (p!0|). 
Now, we rewrite the ZCR (^) in its equivalent (characteristic) form 

fC^ + gC, = VxT (15) 

and decompose the matrix T over the cyclic basis as 

T = PiCu + P2C, + p^VxCu + PiVxC, + P5 V^.C„ + p^VlC, 

+ P7VlCu+PsVlCu, (16) 

where the coefficients pi, . . . ,p8 are functions of a;, t, u, . . . , Ux...x, v, ■ ■ ■ , Vx...x 
and cr. Taking into account the closure equations ([l3| ) and the linear inde- 
pendence of the matrices (|T2|), we obtain from (|T3p and (|TBp the following 
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relations: 



pi = -D^p3 - asps - hpa, p2 = -D^p^ - a^ps - hp&, 
P3 = -Z^^cPs - asPs - Pi = -^xPe - ctePs - Me, 

P5 = -D^p7 - ajps - bjpe, pj = -D^p^ - a^ps - hpQ. 



These relations (|T^), together with (|T3p, determine Pi,P2yP?.-iP4.-iPbyP7y f 
and g in terms of some linear differential operators applied to the two (still 
undetermined) functions p^ and pg- For what follows, we need no expression 
for T: it is sufficient to know that T exists. The expressions for the right-hand 
sides / and g of the represented equations (|ToD are 

h={M + \L + \^K) r, (18) 

where 

A = 9a^ r = f^^^ , (19) 



^/ \P& 



(-6u + v)D^-v^ 



with 



(21) 



Mn = Dl+i -13m + ^t; ) Dl + ( -34m, + ) 



2 y ^ V 2 
3 3 

+ ( -llu^^^ + -^xxx + 33MMa; - -fMx - 3Mfa; ) , (23) 



Mi2 = ^I^^ - uD^ - ^M,, (24) 
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M21 = (-36m, + 8v,) Di 

+ {—l^iUxxx + Sfxxx — 5AuUx + 33vUx + lAuv^ — Gw^) , (25) 



M22 = Dl+(u- ^v^ Dx - ]^Vx. (26) 

If a in (IP is a fixed constant, then ([T8| ) and (p!9D-(p^) give the solution of 
our problem: a continual class of coupled evolution equations, containing two 
arbitrary functions and pg? is represented by (|) with (^. However, a is 
a free parameter in (P), and we need to take into account that the condition 
dh/dX = must be satisfied for (p^S]). Using the series expansion of r 

r = ro + Ari + AV2 + AV3 + AV4 + • ■ ■ (27) 

(if a singularity is suspected in r at A = 0, one may use ( P7D after an in- 



finitesimal shift of A, which does not affect final results), we obtain from ([I8D 
that 



and 



h = Mro (28) 

= Mri + Lro, (29a) 

= Mr2 + Ln + Kro, (29b) 

= Mr3 + Lr2 + is:ri, (29c) 

= Mr 4 + Lr^ + Kr2, (29d) 



This gives us the following implicit solution of our problem: the system of 
coupled evolution equations ([lOD admits the ZCR (^ with the matrix X given 
by (^) if and only if its right-hand side h = {f,g)'^ is determined through 
(|28|) by the set of two-component functions ro, ri, r2, . . . satisfying the three- 
term recurrence ([29|), where the matrix differential operators M, L and K 
are defined by (|2g)-(|2|). 

3 Recursion of systems 

Let us solve the recurrence (PD|), using formal inversion of differential op- 
erators. We say that the operator is formally inverse to a linear dif- 
ferential operator P if P^^P = PP^^ = I, where / is the unit operator. 
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In other words, P~^a denotes any function h of local variables, such that 
a = Pb, if b exists. We have to remind that, in the present case, differential 
operators contain total (not partial) derivatives and act on differential 
functions (in Olver's sense |0, i.e. functions of independent variables, depen- 
dent variables and finite-order derivatives of dependent variables). For this 
reason, D~^0 = (t) with any function 0, whereas {Dl + u) ^0 = and 
{SD^ — AvDx — 2vx) = — |, to give some examples. 



For the operators L and M, defined by (p0|)-([26|), we notice that the 

inverse operators and L"^ exist, but does not exist, and therefore 

we cannot use inversion of K in what follows. We obtain 

ro = -L-^Mri (30) 

from (p9a|) , and 

ri = [KL-^M - L) Mra (31) 

from (129 b|), and notice that the form of (|3Tf) is different from the form of 



(|30|) , because KL ^ for the operators defined by (|20|) - (|26|) . Then we 
obtain 



r2 



(^-K [KL'^M - L) ^M-L^ Mrg (32) 
from ( P9c|) , and find that, fortunately, it is possible to rewrite (|3^) as 

r2 = {KL-^M - L) Mrg (33) 
owing to the condition 

KL-^K = (34) 



satisfied by K (^) and L (|2lD . We see that the form of (p3D coincides with 
the form of (PTj). The same is true, also owing to (|5^), for the relation 
between and r^, which follows from (|29d| ); and, in general, 

ri= {KL-^M ~ Mri+i, i = 1,2,3,.... (35) 

We have to emphasize that the recursion ( ^5)) is caused by the specific prop- 
erty ( p4D of the operators K ( pO] ) and L (pT]) , and that it would be impossible 



to solve the three-term recurrence (29) for generic operators. 



For convenience in what follows, we rewrite (^) as 

ri = M-^LM-^RML-^Mri+i, i = 1,2,3,..., (36) 
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where 

R = MN-^ 



with 

N = ML-^K - L. (; 
The exphcit expressions for and R are 



and 



R ~ ( -^12 I 
R22 J 

with 

= 3Dl - 6m - 3m^D-\ 



i?i2 = [-2Dl + (2m + 3v) Dl + (-4m, + 8m,) DI 
+ (-6m,, + 7m,, + 4m^ - 6mm) L), 
+ (—2m,,, + 2m,,, + 6mm, — 3mm, — 4mm,) 



i?2i = QDl + (6m - 9m) - 3m,L'-\ (^ 

i?22 = [-3Dl + (-18m + 12m) Dl + (-27m, + 18m,) DI 
+ (-21m,, + 14m,, + 12m^ + 12mm - 9m^) 
+ (—6m,,, + 4m,,, + 12mm, + 6mm, + 6mm, — 9mm,) 

+m,D;^,] {3Dl - 4mD, - 2m,) . (^ 
Now, using (^), (|30|) and (p6D, we obtain the following expression for 

h= -R^'ML-^Mrn+u ri = 0,1,2,.... (^ 
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Since the leading-order terms of the operator R (|40|)-(pD are given by 
3DI + --- -iDl + ■ 



and the orders of derivatives of u and v in h must be finite, we can choose 
a sufficiently large n in so that the orders of derivatives of u and v in 
ML~^Mrn+i are less than 2 for this n. In this case, we conclude, taking 
into account the form of M and L (0), that L~^Mr„+i must be a 

function of x and t only, that the orders of derivatives of u and v in Mr^+i 
must be less than 2, and that r^+i = (0, 0)'^. Then we find that 

where the function (t) is arbitrary, and the factor 2 is taken for convenience. 
Finally, (|^) gives us 

h = <j>{t)R-h), \f<j>{t), n = 0,l,2,..., (48) 



or, equivalently, 

/i = i?"Q, n = l,2,3,..., (49) 

because R (0, 0)'^ = (0 (t) Mj,, (t) f^.)'^ with any (t). 

Thus, the problem, formulated in the beginning of Section 2, is solved 
explicitly: the relation (^) determines the right-hand sides h = (/, g)"^ of all 
systems ( pUj ) which admit the linear problem (|^) with the matrix X (^. The 
constructed operator R (|40|)-(P^ turns out to be a recursion operator in the 



Lax sense: it generates a hierarchy of integrable systems which all possess 
Lax pairs with the predetermined spatial part. The system (||) corresponds 
to n = 2 in (^9]) , and the next member of the hierarchy (n = 3) is the 
fifth-order system 

_ 5 40 50 

3 3 3 

10 , 125 

3~^^^^^ ~ ISOMj^Mxx + 4:0VxUxx H -^UxVxx — lOfx^xx 

2 . 5 2 80 2 20 

-|- 120m U^; — A\)uVUx H 1^ U Vx H tit'Wa;, 

3 3 3 

5 55 

- ...... - m^UxUxx + 80..^.. + mUxVxx - ^VxVxx 

35 

+ IQQuvUx — 'iOv'^Ux + 'iOu^Vx — GOuvVx + —v^Vx- (50) 
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4 Recursion of symmetries 



The operator R, defined by (^)-(|4^), is also a recursion operator with re- 
spect to symmetries of the coupled KdV equations @. In order to see this, 
we have to prove that R satisfies the following condition [|]: 

DtR=[H,R], (51) 

where H is the Frechet derivative of the right-hand side of the system 

with 

Hu = 4:Dl + {-12u + v) + (-12m, + 2v,) , 
Hi2 = -Dl + 2mD, + M„ 
H21 = 9Dl - l2vD, - Qv,, 

H22 = -2Dl + (-6m + Av) D, + (-12m, + 4m,) . (53) 
It is convenient to study the condition (|5TD in its equivalent form: 



HM - DtM = R {HN - DtN) , (54) 

which follows from (^) through the definition of R (|37| ) and the identity 
DtN^^ = —N'^ {DtN)N~^. Direct computations show that the condition 
(M) is satisfied identically by H m)-(m), M (|22|)-(|2B|), R (m-(M), N 



and the expressions (||) for Ut and Vt- Since the condition (|51|) is satisfied, the 
system of coupled KdV equations (|^) possesses infinitely many symmetries: 
the recursion operator R generates a 'new' symmetry (but not necessarily a 
local one) from any 'old' symmetry. 



5 Conclusion 

The recursion operator R (^0|)-(^), constructed in this paper for the new 
integrable system of coupled KdV equations (||), is characterized by unusual 
structure of its nonlocal part, which contains not only the conventional oper- 
ator D^^ but also the strange nonlocal operator {3Dl — 4mD, — 2m,) \ Such 
a phenomenon has not been observed in the literature as yet. This opera- 
tor {3D^ — 4:vDx — 2m,) ^ cannot be expressed in terms of a finite number 
of operators D~^, if one uses only local variables, as we did throughout the 
paper. However, if we introduce the nonlocal variable w : v = 3w,,/w, we 
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can represent the obtained recursion operator in a more conventional form, 
using the factorization {SD^ — 'ivD^ — 2f^.)~^ = ^w"^ D^^w^'^ D~^w~'^ D^^w'^ . 

We found the recursion operator of the system by the technique of 
cychc bases of ZCRs. It is interesting if this recursion operator can be ob- 
tained by different methods as well, e.g. by those developed in |0, ||, |^. 

We have seen that the relation (|49|) generates local expressions for h at 
n = 1, 2, 3. It is easy to check that the seventh-order system of this hierarchy 
(n = 4) is also a local expression. Nevertheless, we have not proven as yet 
that this is the case for any n. 

Though our method produces the recursion operator in the quotient form 
MN~'^ (13), the operators M (|2D-(|6D and (H) are not Hamiltonian 
operators. It is still unknown if the obtained recursion operator is hereditary 
and if the system (H) admits a bi-Hamiltonian structure. 
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